Dimensional crossover and the magnetic transition in electron doped manganites 
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We introduce a microscopic model for electron doped manganites that explains the mechanism of 
the observed transition from G-type antiferromagnetic (G-AF) to G-type antiferromagnetic (G-AF) 
order under increasing doping by double exchange mechanism. The model unravels the crucial role 
played by e 9 orbital degrees of freedom and explains the observed metal-to-insulator transition by 
a dimensional crossover at the magnetic phase transition. The specific heat and the spin canting 
angle found for the G-AF phase agree with the experimental findings. As a surprising outcome 
of the theory we find that spin canting is suppressed in the G-AF phase, in agreement with the 
experiment, due to the Fermi surface topology. 
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I. INTRODUCTION 

Doped transition metal compounds cxibit a num- 
ber of very interesting and fascinating phenom- 
ena, including high temperature superconductivity in 
cupratcsA or pnictides^ and colossal magnetoresistance 
in manganites^ The hole-doped perovskite manganites 
are well known for their complexity, originating from 
the competing tendencies in their rich phase diagrams^ 
with a variety of spin, charge and orbital orders. Several 
phases are energetically close to each other and small 
changes of electronic parameters or electron concentra- 
tion result in phase transitions^^ Hole doping, realized 
e.g. in Lai_j,(Ca,Sr)j,Mn03 for increasing y, triggers 
then a transition from an antiferromagnetic (AF) to a 
ferromagnetic (FM) phase by the double-exchange (DE) 
mechanism^ that became a crucial idea to explain the 
onset of the metallic FM phase and electronic transport. 

The DE model is one of the most widespread mod- 
els of ferromagnetism. It describes the kinetic energy 
of electronic charge carriers coupled by Hund's exchange 
to localized ti g spins interacting by AF superexchange. 
A pioneering work by dc Gcnncs^ suggested that the 
competition between FM DE interaction that arises in 
doped systems and AF superexchange between ti g spins 
leads to spin-canted AF order. Extensive studies of these 
ideas led to progress in the theoretical understanding 
of the phase diagrams of hole doped manganites 
while electron doped systems were not analyzed care- 
fully until now. In recent years the DE model was used 
to study, inter alia, striped multifcrroic phases in doped 
manganites ) 10 ' 11 evolution of spin, orbital and charge or- 
der at interfaces^ valence fluctuations in magnetite jA 3 -, 
novel spin phases on frustrated lattices^ and unuasual 
spin superstructure observed in several vanadium spinels 
with degenerate t^g orbitals*^ In all these cases e g or- 
bital degeneracy plays a crucial role and might lead to 
anisotropic AF phases: either C-type AF (C-AF) or A- 
type AF (A-AF) phased depending on the electronic 
filling and the strength of AF superexchange^ 



Usually, however, accurate studies of the DE model 
were hindered by the presence of strong intraorbital 
Coulomb interaction U in a lattice distorted by strong co- 
operative Jahn- Teller effecti^^ A qualitatively different 
situation arises in an ideal cubic perovskite SrMn03, with 
e g bands being completely empty. Here the only interac- 
tion is AF superexchange between 5 = 3/2 spins at Mn 4+ 
ions. They are generated by t2 g charge excitations 1 ^ 
and stabilize the isotropic three-dimensional (3D) G-type 
AF (G-AF) phased Weakly doped SrMn0 3 represents 
an ideal situation for the DE model as the e g electrons 
hardly interact with each other at low doping x and 
strong on-site Coulomb interactions between them may 
be neglected. In this case the DE model can be almost 
exactly solved under a single approximation, i.e., neglect- 
ing weak quantum fluctuations of large 5 = 3/2 spins. 

Early studies of electron doped Ca 1 _ a ,/ 2 Ce a: /2Mn03 
systems gave a magnetic transition to the G-AF phased 
The magnetic order in doped CaMn03 was analyzed 
within the density-functional theory calculations that 
give instead a transition from the G-AF to A-AF 
phased Recently single crystals of Si'i_ T La x Mn03 and 
Sr 1 _ ;E /2Ce z: /2Mn03 were synthesized and investigated^ 
showing that the transition to the G-AF phase occurs 
in both systems, in agreement with the prediction by 
Maezono et al.^ and the magnetic phase diagram is re- 
markably universal as a function of doping x. 

In electron doped manganites chemical doping of Sr 2+ 
with La 3+ (Ce 4+ ) ions generates e g electronic carriers 
at Mn 3+ ions. Following the standard picture of de 
Gennes;£ one expects spin canting on AF bonds. Until 
recently only one exception from this rule was reported 
— spin canting is absent in the A-AF phase when elec- 
trons occupy x 2 — y 2 orbitals and the hopping along AF 
bonds between FM ab planes is blocked^ In this context 
the recent experimental result of no canting in the G-AF 
phase 2 ^ is puzzling and requires theoretical explanation, 
as hopping on the AF bonds in ab planes is there finite 
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Hole doped Lai-ySr^MnOa manganitcs with large e g 
electron density remain insulating up to y ~ 0.17^ as the 
metallic behavior is hindered here by formation of orbital 
polaronsi 23 i 24 . This is in sharp contrast with the elec- 
tron doped systems with the G-AF metallic phase found 
in Sri-aXa^MnOa by Sakai et alM. already at a very 
low electron doping x = 0.01. Surprisingly, when dop- 
ing increases further to x ~ 0.04 this 3D weakly metal- 
lic phase changes to an insulating phase with quasi one- 
dimensional (ID) orbital order of partly occupied 3z 2 — r 2 
orbitals, accompanied by tetragonal lattice distortion and 
the C-AF spin order. 

In this paper we investigate the DE model for doped 
e g electrons, extended here by the coupling to the tetrag- 
onal lattice distortion which occurs in the C-AF phase. 
We derive the exact electronic structure, show a remark- 
able difference between the Fermi surface topology in the 
above magnetic phases, and demonstrate that it explains 
the absence of spin canting in the C-AF phase. The cal- 
culated values of the specific heat coefficient 7 and the 
critical doping x c for the G-AF/C-AF phase transition 
are in quantitative agreement with experiment. 

The paper is organized as follows. In Sec. [TT| we in- 
troduce the microscopic model and specify its parame- 
ters. The model is solved in Sec. IIIIl where we analyze 
the remarkable difference between the electronic struc- 
ture found in two magnetic phases. This leads to the 
magnetic transition which occurs under increasing dop- 
ing as we show in Sec. IIVI The paper is summarized in 
Sec. El 

II. THE DOUBLE EXCHANGE MODEL 

Wc consider the DE model for degenerate e g 
electrons,— extended by the coupling to the lattice, 

ij,atf3 : <T i 

+ JJ2 Si ■ Sj - guJ2(n lz ~ n lx ) + ^NKu 2 . (1) 
(ij) » 

Here the c\ aa operator creates an electron with spin a =t 

, I in orbital a = z, x at site i, and rii a = J2 a c iau c jaa m 
electron density in orbital state a at site i, and N is the 
number of sites. Two e g orbitals are labeled as 

z = (3z 2 - r 2 )/V6, x = (x 2 - y 2 )/V2, (2) 

and correspond to two components of the pseudospin r = 
1/2. The hopping elements between a pair of nearest 
neighbors (ij) along the axis 7 = a, b 1 c are: 

where t in an effective (dda) hopping element between 
two 3z 2 — r 2 orbitals along the c cubic axis. Following the 



estimates of t performed within the polaronic picture 2 ^ 
and using the optical spectroscopy for LaMnOaj 2 ^ we fix 
the effective (dda) hopping element at t = 0.4 eV. 

Hund's exchange between e g and t2 g electrons at Mn 4+ 
ions Jh — 0.74 eV is here the largest parameter. It 
aligns the spin of an e g electron with the localized 
tig spin Si at site i. In this way the electronic struc- 
ture depends crucially on the magnetic order of localized 
t2g spins {Si} that interact by the nearest neighbor AF 
superexchange J > 0. It follows from virtual charge ex- 
citations (if dj ^ dfd 2 along nearest-neighbor bonds (ij) 
in the presence of strong on-site Coulomb interaction U 
between t-2 g electrons.— 

The tetragonal distortion u is finite only in the C-AF 
phased Here we define it as proportional to a difference 
between two lattice constants a and c along the respective 
axis, u = 2(c — a)/(c + a). Increasing lattice distortion 
(u > 0) causes increasing e g orbital splitting AEjt de- 
termined by the Jahn- Teller energy, Ejt = g 2 jlK ~ 0.2 
eVi 2 ^ Furthermore, using the experimental data of Rcf. 
|22| . we arrived at a semiempirical relation u ~ 3a;/20. 



III. ELECTRONIC STRUCTURE 

Following dc Gcnnes^ we determined the electronic 
structure for the G-AF phase by considering spin canting 
by angle 6 caused by electron doping. Hund's exchange 
at each site favors aligned ti g and e g electron spin states. 
To solve the electronic structure we introduce the fermion 
{flaa} operators in a canted spin background using the 
local coordinate frame, for instance, 

4^ = co S (e/2)fl t -sm(e/2)f^. (4) 

The canting of each spin by angle 9 £ [0, n/2] (other 
states are equivalent) corresponds for 9 = (tt/2) to the 
AF (FM) spin order. The canting modulates the hopping 
elements between neighboring sites, and one finds an an- 
alytic solution for two (majority) e g bands in the G-AF 
phase for spins being parallel to the t2 g spins at each site 
due to Hund's exchange, 

4 = -^/(tsin0 7 ± + J H S) 2 + (icos0 7 ±) 2 , (5) 

where 

7k = ^2 cos k x ± I ^2 cos2 fc A 

A A 

1 1/2 

— cos k x cos k y — cos k y cos k z — cos k z cos k x > , 

(6) 

and the summations in Eq. (|6]) are over three cubic direc- 
tions in momentum space, i.e., A = x,y,z. The remain- 
ing e g bands correspond to excited states of e g electrons 
with spins antiparallcl to those of localized t^ g electrons 
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FIG. 1: (Color online) Band structure in the G-AF phase 
along the high symmetry directions. Spin majority (minority) 
bands are shown by solid (dashed) lines. Parameters: t — 0.4 
eV, Jh — 0.74 eV, x — 0. Inset shows the Fermi surface at 
low doping (heavy lines along V — X , R — M, and equivalent 
directions) and the special points: r = (0, 0, 0), X — (tt, 0, 0), 
M = (tt, tt, 0), R = (tt, tt, tt), Z = (0, 0, tt). 



— they are given by — e k and separated from the lower 
bands (J5]) by a large gap of 2JhS at 7 k =0. 

In the undoped G-AF phase of SrMn03 spin canting 
is absent (9 = 0) and the electronic structure consists of 
four Slater bands ±e k , see Fig. [TJ Both e g orbitals are 
equivalent here and equally contribute to each electronic 
band. The lowest energy is obtained along the Y — X 
and equivalent directions in the Brillouin zone and here 
electrons are doped. By considering the electron energy 
close to these lines one finds that the Fermi surface in 
the regime of low electron doping, shown in the inset in 
Fig. [TJ consists of narrow cylinders along equivalent main 
directions in the reciprocal space. 

The canting angle in the G-AF phase is given by 



sint 



tJ H S 



V(3i) 2 + ( J H Sf 12J5 2 



(7) 



and increases linearly with x in the low doping regime. 
Using an empirical formula for the Neel temperature of 
a Heisenberg antiferromagnet^ we obtained J = 3.79 
meV for SrMn0 3 from the experimental T N ~ 230 KM 
For the parameters of Fig. [TJ one finds from Eq. ([7]) that 
the canting angle is rather small, sin (9 ~ 2.7x, and quite 
close to the experimental value, 8 exp ~ 2xM 

The electronic structure {ek} for the G-AF phase is 
quite distinct from the one in the G-AF phase and was 
determined numerically. It depends on doping x as the 
tetragonal lattice distortion increases with it. We have 
verified that this dependence is weak for the realistic pa- 
rameters, i.e., with g ~ 3.0 eV and K ~ 20 eVM A rep- 
resentative result is shown in Fig. [2]for x = 0.05. As the 
FM order along the c axis breaks the cubic symmetry, the 
orbitals decouple from each other in the electronic struc- 
ture and form two independent subsets of bands. This is 
best visible along the Y — Z directions where one band 



FIG. 2: (Color online) Band structure in the C-AF phase 
along the high symmetry directions. Spin majority (minority) 
bands are shown by solid (dashed) lines. Parameters: t — 0.4 
eV, Jh = 0.74 eV, g = 3 eV, x = 0.05. The special points 
as in Fig. [TJ Inset shows the Fermi surface pockets at low 
doping near the X and Y points. 



has a large dispersion close to At, and the other one is 
dispersionless (Fig. [2J. In addition, the former band has 
also a very weak dispersion along two (equivalent) Y — X 
and Y — Y directions, with two minima at the X and Y 
points that arise due to the weak mixing between z and 
x orbitals by hopping in ab planes, but this hopping is 
almost fully quenched due to the AF order. 

Surprisingly, we found no spin canting in the G-AF 
phase. Here the spins form an ideal staggered AF struc- 
ture in the ab planes, while they are aligned in FM chains 
along the c axis. The reason of this behavior is the elec- 
tronic structure with the lowest energy at the X and Y 
points and the Fermi surface consisting of two pockets 
that are centered around these two equivalent points at 
low doping, see the inset in Fig. [5] For this Fermi sur- 
face topology no kinetic energy can be gained by spin 
canting, since the resulting nearest-neighbor hopping be- 
tween z orbitals in ab planes, jts'm8(cosk x + cosk y ), 
vanishes at the X and Y points. Thus, the observed sup- 
pression of spin canting in the G-AF phase follows from 
the abrupt change of the Fermi surface topology at the 
magnetic transition. 

The consequences of flavor separated band structure 
in the G-AF phase are remarkable. First of all, electrons 
doped in SrMnOs have pure z orbital character. Second, 
this implies that they can propagate only along the c 
axis where localized spins are FM, but their motion is 
blocked in the ab planes by the AF order, in agreement 
with DE. Altogether this leads to a dimensional reduction 
of electronic transport to ID chains along the c axis. 

A qualitatively different character of the electronic 
structure in both AF phases is visible in the orbital- 
resolved density of states (DOS) N z (ui) and N x (ui) in the 
spin majority bands, see Fig. [3l In the G-AF phase both 
orbital DOSs arc the same and have low bandwidth of 
~ 0.5 eV due to the AF order in all three cubic directions 
[Fig. [3{a)]. At the lower band edge one finds a character- 
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FIG. 3: (Color online) Densities of states for e g electrons, 
N z (co) and N x (oj), in the spin majority bands as obtained for: 
(a) the G-AF phase (N z (oj) and N x (u)) are here degenerate) 
with a singularity arising from the AF order at uj = —1.11 eV, 
and (b) the C-AF phase; here the density of z, N z (lj), and 
x, N x {ui), electron states is shown by solid (red) and dashed 
(blue) line. Parameters: t = 0.4 eV, J H = 0.74 eV; for C-AF 
phase in addition g — 3 eV and x = 0.05. 



istic step which reflects the quasi-two-dimensional DOS 
of e g electrons^! A large gap between the spin majority 
and spin minority states in the AF phase generates a van 
Hove singularity at the upper edge of the DOSs. 

From the DOS in the G-AF phase, see Fig. [3ja), we 
have estimated the coefficient 7 ~ 5 mJ/mol-K 2 in the 
specific heat cy ~ 7T at x = 0.01 which is indeed very 
close to the experimental value 7 eX p — 5.4 mJ/mol-K 2 , 
and weakly increases with doping x as in experiment^ 
Thus we find that almost no polaronic renormalization 
of the electron mass occurs in Sr 1 _ a .y 2 Ce a; /2Mn03 and 
Sri^La^MnOs. 

Unlike in the G-AF phase, the orbital flavors separate 
in the G-AF phase into a wide quasi-lD z-subband of 
width ~ At and a narrow (empty) x-subband, similar to 
the majority spin bands in the G-AF phase. The former 
subband gives a characteristic ID DOS N z (u>), with two 
maxima at both band edges. Due to the AF order in the 
ab plane, the separation of the e g orbital flavors is almost 
perfect here. Large DOS N z (u) at the edge of the quasi- 
lD band [Fig. |3[b)] is here sufficient for self-trapping of 
electrons^ within spin polarons, as is doped CaMn03j22 
In the present case we believe that the polarons arc quasi- 
lD, in agreement with the nature of underlying G-type 
structure and showing a similar behavior to orbital po- 
larons in a ID model of hole-doped manganites<2i We 
thus expect a strongly anisotropic optical conductivity 
in the insulating G-AF phase: a a b(uj) with a large gap, 
and a pseudometallic a c (uj) with small charge gap. 



IV. MAGNETIC PHASE TRANSITION 

Now we show that the DE triggers the magnetic transi- 
tion in an electron doped system by comparing the total 
energies in both magnetic phases. In the undoped sys- 



0.00 




X 



FIG. 4: (Color online) Total energies of the G-AF (Egaf) 
and G-AF (Ec af) phase as obtained for increasing electron 
doping x. Parameters: t = 0.4 eV, Jh = 0.74 eV, g — 3 eV, 
K — 20 eV, and J = 3.79 meV. Inset shows the phase diagram 
in the (t,x c ) plane obtained for: g = (circles, dashed line), 
and g = 3 eV (diamonds, solid line). 



tern [x — 0) the ground state is the G-AF structure, with 
low magnetic energy of Eq AF = —3JS 2 per site, and this 
energy increases to — 3JS 2 cos 28 for finite doping x > 0. 
In the collinear G-AF phase this energy is enhanced to 
Eq AF = —JS 2 , as the superexchange bonds along the 
c axis are here frustrated. Simultaneously, however, in- 
creasing doping gives a different kinetic energy in the 
above phases, obtained by summing up the energies of 
occupied states for a given doping x. The lower edge of 
the DOS is much lower in case of the G-AF phase than 
for the G-AF phase; for instance at x = 0.05 these min- 
ima are at luq ~ —1.63 eV and wc — —2.07 eV, see Fig. 
O One finds the energy gain due to full electron hopping 
~ t along the c axis in the ID electronic structure of the 
G-AF phase, while only a small fraction of t contributes 
to the kinetic energy in the isotropic 3D G-AF phase. 
Therefore, the total energy of the G-AF phase decreases 
faster with increasing doping x, see Fig. 2J and finally 
this phase becomes more stable than the G-AF one for 
x > x c , as observed in the experiment^ We note that 
small spin canting in the G-AF phase has almost no effect 
on the phase boundary as the gain in the kinetic energy is 
almost cancelled by the loss of the superexchange energy. 

The crucial role played by the DE mechanism is visible 
in the phase diagram of electron doped manganites shown 
in the inset of Fig. [4] For increasing hopping t both 
phases have increased kinetic energy gains, but a larger 
energy gain is found in the G-AF phase. As a result, the 
critical concentration x c is lowered when t increases, but 
one finds that x c lies still in a relatively narrow window 
of doping 0.03 < x c < 0.04 in a broad range of realistic 
values of t e (0.3, 0.6) eV. 

The magnetic transition is first order and occurs at a 
somewhat lower doping x c when the tetragonal distortion 
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which stabilizes the C-AF phase is considered, see the in- 
set of Fig. 01 However, this energy gain is rather small 
compared with the kinetic energy difference between the 
two magnetic phases, so we conclude that the tetragonal 
distortion has only a weak stabilizing effect for the ID 
3z 2 — r 2 structure in the C-AF phase, while the transition 
itself follows from the DE interaction. We remark that 
this scenario does not apply to La x Cai_ ;r Mn03 nanopar- 
ticles where phase separation into AF and FM domains 
takes placed 

V. SUMMARY 

Summarizing, we have introduced the microscopic 
model that explains the mechanism of the magnetic tran- 
sition in electron doped manganites from canted G-AF 
to collinear C-AF phase at low doping x ~ 0.04. We 
demonstrated that the DE supported by the cooperative 
Jahn- Teller effect leads then to dimensional reduction 



from an isotropic 3D G-AF phase to a quasi-lD order 
of partly occupied 3z 2 — r 2 orbitals in the C-AF phase. 
This prediction of the theory can be verified by future 
angle resolved photoemission experiments as the shape 
of the Fermi surface in the G-AF and C-AF phases is 
radically different. 

The presented theory explains as well the absence of 
spin canting in the C-AF phase by the Fermi surface 
topology. This is a subtle effect as electron doping oc- 
curs here near the X and Y points in the Brillouin zone, 
and the kinetic energy cannot be gained in spin canted 
structure. 



Acknowledgments 

We thank Lou-Fc' Fciner and Peter Horsch for insight- 
ful discussions. A.M.O. acknowledges support by the 
Foundation for Polish Science (FNP) and by the Polish 
Ministry of Science under Project N202 069639. 



1 P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006); M. Vojta, Adv. Phys. 58, 699 (2009). 

2 J. Paglione and R. L. Greene, Nat. Phys. 6, 645 (2010). 

3 E. Dagotto, T. Hotta, and A. Moreo, Phys. Rep. 344, 1 
(2001); E. Dagotto, New J. Phys. 7, 67 (2005). 

4 Y. Tokura, Rep. Prog. Phys. 69, 797 (2006). 

5 A. Weifie and H. Fehske, New J. Phys. 6, 158 (2004). 

6 C. Zener, Phys. Rev. 82, 403 (1951). 

7 P. G. de Gennes, Phys. Rev. 118, 141 (1960). 

8 R. Maezono, S. Ishihara, and N. Nagaosa, Phys. Rev. B 
58, 11583 (1998). 

9 T. Hotta, A. L. Malvezzi, and E. Dagotto, Phys. Rev. B 
62, 9432 (2000); A. M. Oles and L. F. Feiner, ibid. 65, 
052414 (2002); D. V. Efremov, J. van den Brink, and D. I. 
Khomskii, Nature Mat. 3, 853 (2004). 

10 S. Dong, R. Yu, J.-M. Liu, and E. Dagotto, Phys. Rev. 
Lett. 103, 107204 (2009). 

11 S. Liang, M. Daghofer, S. Dong, C. §en, and E. Dagotto, 
Phys. Rev. B 84, 024408 (2011). 

12 G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 101, 216804 
(2008). 

13 R. J. McQueeney, M. Yethiraj, S. Chang, W. Montfrooij, 
T. G. Perring, J. M. Honig, and P. Metcalf, Phys. Rev. 
Lett. 99, 246401 (2007); 100, 069901(E) (2008). 

14 Y. Motome and N. Furukawa, Phys. Rev. Lett. 104, 106407 
(2010); S. Kumar and J. van den Brink, ibid. 105, 216405 
(2010). 

15 G.-W. Chern and C. D. Batista, Phys. Rev. Lett. 107, 
186403 (2011). 

16 In the C-AF (A-AF) phase the spin order parameter is 
staggered in ab planes (along the c axis), while in the G- 
AF phase it is staggered in all three cubic directions. 

17 J. van den Brink and D. Khomskii, Phys. Rev. Lett. 82, 



1016 (1999). 

18 I. Leonov, D. Korotin, N. Binggeli, V. I. Anisimov, and D. 
Vollhardt, Phys. Rev. B 81, 075109 (2010). 

19 L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 

20 Z. Zeng, M. Greenblatt, and M. Croft, Phys. Rev. B 63, 
224410 (2001); E. N. Caspi, M. Avdeev, S. Short, J. D. Jor- 
gensen, M. V. Lobanov, Z. Zeng, M. Greenblatt, P. Thiya- 
garajan, C. E. Botez, and P. W. Stephens, ibid. 69, 104402 

(2004) . 

21 H. Tsukahara, S. Ishibashi, and K. Terakura, Phys. Rev. 
B 81, 214108 (2010). 

22 H. Sakai, S. Ishiwata, D. Okuyama, A. Nakao, H. Nakao, 
Y. Murakami, Y. Taguchi, and Y. Tokura, Phys. Rev. B 
82, 180409 (2010). 

23 R. Kilian and G. Khaliullin, Phys. Rev. B 60, 13458 
(1999). 

24 M. Daghofer, A. M. Oles, and W. von der Linden, Phys. 
Rev. B 70, 184430 (2004). 

25 A. M. Oles, P. Horsch, G. Khaliullin, and L. F. Feiner, 
Phys. Rev. B 72, 214431 (2005); N. N. Kovaleva, A. M. 
Oles, A. M. Balbashov, A. Maljuk, D. N. Argyriou, G. 
Khaliullin, and B. Keimer, ibid. 81, 235130 (2010). 

26 G. S. Rushbrooke and P. J. Wood, Mol. Phys. 1, 257 
(1958); M. Fleck, M. G. Zacher, A. I. Lichtenstein, W. 
Hanke, and A. M. Oles, Eur. Phys. J. B 37, 439 (2004). 

27 L. F. Feiner and A. M. Oles, Phys. Rev. B 71, 144422 

(2005) . 

28 G. Wellein and H. Fehske, Phys. Rev. B 58, 6208 (1998). 

29 H. Meskine, T. Saha-Dasgupta, and S. Satpathy, Phys. 
Rev. Lett. 92, 056401 (2004). 

30 Y. Wang and H. J. Fan, Phys. Rev. B 83, 224409 (2011). 



